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Abstract 

^D , We present the results of a statistical analysis of the Doppler shifts and the asymmetry 

' parameters of V profiles of the Fe I 630.25 nm line produced by 2D MHD simulations 



of solar granulation. The realism of the simulations tested using the magnetic ratio of 
Fe I 524.71 and 525.02 nm lines. The Stokes spectra were synthesized in snapshots with 
a mixed polarity field having a mean magnetic fiux density of 0.2 niT and mean unsigned 
"^ ■ field strength of 35 mT. We found that downflows with a velocity of 0.5 km s~^ predomi- 

nate, on the average, in areas with some network magnetic elements at the disk center. In 
H ! separate strong fluxtubes average velocity is equal to 3 km s~^ and the maximum velocity 

^ ' is 9 km s~^. In weak diffuse magnetic fields upfiows dominate. Their average velocity 

Cd . is 0.5 km s"^ and maximal one is 3 km s~^. The ^-profile asymmetry depends on the 

spatial resolution. The V profiles synthesized with high spatial resolution (35 km) have 
average amplitude and area asymmetries —1%, 1%, respectively. The asymmetry scatter 
is ±70% for weak profiles and ±10% for strong ones. The profiles with low spatial resolu- 



t~^ , tion (700 km) have average amplitude and area asymmetries 3%, -2%, respectively. Low 

^^ I spatial resolution is a reason why the amplitude asymmetry is always positive and greater 

than the area asymmetry in V profiles observed. We found weak correlation between the 



asymmetry of V profiles and velocity. Upfiows cause negative asymmetry, on the aver- 

f^ I age, and downfiows cause positive asymmetry. We examined center-to-limb variations of 

vertical velocity in magnetic elements. Beginning from cos 9 = 0.9, the average velocity 

abruptly increases from 0.5 to 2 km s~^ and then slightly varies closer to the limb. We 

k><( \ found nonlinear oscillations of vertical velocity with power peaks in the 5-minute and 

^ ' 3- minute bands. This nonlinearity is caused by magnetic field strength fiuctuations in 



fiuxtubes. The Doppler shifts and asymmetry parameters obtained for space- average V 
profiles are consistent with results of FTS observations as well as with other observations 
made with higher spatial resolution. 

1 Introduction 

Magnetic fields on the Sun exhibit fine structure which is conditioned by the fine structure 
of the solar surface in general and its active regions in particular [511]. The first supporting 
evidence for the fine structure of magnetic fields was found in the 1960s, initially in strong 
fields in spots and later in quiet regions [30]. Magnetic fields outside spots were found to be 



unevenly distributed over the surface, in the form of dense magnetic fluxes in very small areas 
(fa 1"). Kilogauss strengths of these compact magnetic features were initially measured in 
faculae and in the supergranular network by the line-ratio method flEl SH US]- These results 
were confirmed more than once in later studies. According to [H], the typical size of magnetic 
concentrations is 100-300 km, typical strength is 100-200 mT, the velocity does not exceed 
1.2 km s~^ and is 0.5 km s~^ on the average. The field strength scatter is very small. The 
authors of [T6] concluded that 90 percent of the photospheric magnetic field is concentrated 
in compact kilogauss structures and 10 percent is found in quiet regions with a mean flux 
density of 0.2-0.3 mT. Further investigations revealed that fine compact magnetic features 
have a filamentary structure and look like tubes with their diameter growing with height in the 
atmosphere (the mushroom or canopy effect). They are tilted most probably at an angle of 10° 
pT| 121] . The temperature in them is higher in the middle and upper photosphere and is lower 
in the lower photosphere as compared to the quiet Sun, and the chromospheric temperature 
growth begins in deeper layers as compared to the quiet Sun [38j. Analyses of all available 
data revealed a wonderful property of these magnetic structures — despite some distinctions 
in size, field strength, and orientation, they all are much alike in form and in their magnetic 
and temperature properties. Because of this, they were pooled to form a separate class of solar 
magnetic structures and were called the small-scale magnetic elements or magnetic flux tubes 
[2^ [55] . The number of flux tubes on the solar surface increases with decreasing tube size. 
Although each magnetic element adds very little to the general solar radiation flux, millions 
of magnetic elements not only compensate for the energy blocked in spots but even produce 
some additional radiation (about 0.1 percent |39]). Magnetic flux tubes can form clusters of 
various sizes. These clusters are assumed to form larger magnetic structures such as the active 
supergranular network, faculae, and plages. At present the smallest (120 km) magnetic elements 
resolved in observations are the so-called bright points of magnetic network [3] . 

Concurrent with the refinement of observation techniques, theoretical models of magneto- 
convection in quiet regions on the Sun have been developed. Much research was devoted to 
numerical simulation of magnetic elements, their formation and interaction with convective 
motions (see review [2S])- 

In this study we use two-dimensional MHD models developed by Gadun [HI IHl 125] • Unlike 
other researchers, he used an original approach to the solution of the equations of radiation 
magnetohydrodynamics [7] which describe a compressible, gravitationally stratified turbulent 
medium. The magnetic field is described by the vector potential, so that the divergence of 
magnetic field strength is always zero in the simulation region. The initial magnetic field is 
specified by a bipolar magnetic configuration. The first sequence of 2D MHD models [B] of com- 
pletely nonstationary solar magnetogranulation extended over a long time period (two hours of 
solar time). With such models, the evolution of solar magnetic elements, their structure and 
dynamics could be studied in detail on scales much smaller than the spatial resolution of the 
present-day observations. Three-dimensional models (e.g., [H]) are certain to give more realistic 
flow patterns, while two-dimensional models more adequately represent small-scale phenomena, 
and they still remain useful in studying the properties of magnetic elements. Analyses of 2D 
simulations of convection [H [2Z] and magnetoconvection [21 [21 [131 US] suggest that 2D models 
reproduce quite well many features of 3-D convection. The major results of 2D magnetoconvec- 
tion simulations on granulation scales obtained by Gadun were presented in |6j. Investigations 
into the mechanism of formation and destruction of flux tubes and into their stability condi- 
tions revealed that thermal flows are the most important factor in the evolution of small-scale 
magnetic fields. Fragmentation of large-scale thermal flows can result in the formation of 
vertical magnetic elements from horizontal magnetic surface flows. This previously unknown 
mechanism of flux tube formation was called the surface mechanism [71 [9]. The mechanism of 



formation of kilogauss flux tubes through convection collapse was also investigated, and some 
spectropolarimetric manifestations of this process were found [SS]- MHD models [B] pointed 
to the emergence, reconnection, and recycling of magnetic flux near the surface, and these 
processes were demonstrated in [26]. Nonstationary MHD models p] can serve as an analog 
of some observed photospheric regions with average unsigned strength of magnetic fields of 
40-50 mT. These heterogeneous models, which are more realistic than two-component models, 
were also used to calculate the Stokes profiles of photospheric lines [321 [Ml [26] . In particular, 
the analysis of the V profiles of the IR line Fe I A 1564.8 nm [M] gave a magnetic field intensity 
distribution which is in good agreement with observations of magnetic fields in quiet regions 
[19] and which confirms the two-component structure of intra-network fields. 

The magnetic field in the solar plasma is closely associated with the velocity field, and 
this association inside sunspots as well as in the quiet solar atmosphere still is one of the 
most difficult problems in solar physics. The freezing of magnetic field in the moving plasma 
(or its rigid connection with velocity field) can be weakened due to the fine structure of this 
magnetic field [30]. That is why the study of dynamic processes in magnetic elements, and 
inside strong fiux tubes in particular, by the magnetoconvection simulation techniques is of 
current interest, inasmuch as the magnetic elements still are not spatially resolved because of 
their small horizontal size. Numerical MHD simulations P| 113] suggested that a large variety of 
dynamic processes should be expected on scales of several hundred kilometers inside and outside 
fiux tubes. Recent highly accurate polarimetric observations in the line Fe I A 630.25 nm with 
a spatial resolution better than 700 km [36] yielded some statistical relations between the 
observed Stokes V^-profile parameters which are indicative of substantial systematic motions 
in the magnetic elements in the network, faculae, and quiet regions. We believe that the 
reproduction of some observed spectropolarimetric effects [201 IMl ES] with the use of 2D MHD 
models [6] is of obvious interest. 

The main purpose of this study is to determine the dynamical characteristics of magnetic 
elements in the solar photosphere. Below we give some principal results of observations of the 
magnetic element dynamics, a description of the 2D MHD models we use, and the results of 
the application of the line-ratio method to these models. We also analyze our calculations of 
the Stokes V profiles, their shifts, asymmetries, and center-to-limb variations on scales smaller 
than or comparable to the spatial resolution of present-day observations. 

2 Observed shifts and asymmetries of V profiles 

Studies based on complex magnetographic observations of magnetic fields and line-of-sight 
velocities revealed that the vertical velocities are equal to zero, on the average, at the chromo- 
spheric network boundaries as well as inside network cells [50], while the line-of-sight velocities 
determined from spectropolarimetric observations were contradictory [22] . Early extensive po- 
larimetric observations with Fourier spectrometers (FTS) with high spectral resolution and low 
spatial and time resolutions (about 10" and 30 min) detected no significant systematic motions 
(above 0.25 km s~^) inside fiux tubes [IQ]. This result differed essentially from a mean velocity 
of 0.5 km s~^ derived earlier in [H]. Only observations with high spatial resolution (^ 1") 
32 detected a small predominance of downflows with a mean velocity of 0.2 km s~^ and the 
dependence of this velocity on the size of magnetic elements. In regions with small filling factor, 
downfiows were stronger and the velocity scatter increased from —0.5 km s~^ to 1.5 km s~^. 
The mean velocity measured in [11] was about 0.8 km s~^, and the scatter increased with de- 
creasing amount of polarization. Theoretical simulations of the interaction between convection 
and magnetic field [21 [13] also confirmed the existence of systematic motions outside and inside 



flux tubes on small spatial and temporal scales. A most impressive evidence that dynamic 
processes are stronger inside magnetic elements was given in recent studies [121 ES] • Systematic 
downward motions observed with a spatial resolution of 0.8-1" were found to have mean ve- 
locities of 0.5-0.7 km s~^. In magnetic elements with small filling factor the velocities were as 
high as ±5 km s~^, while they were smaller in elements with large filling factor or in clusters. 
Some evidence for large horizontal velocities of about 2 km s~^ in magnetic elements was given 
in [51]. 

In addition to the zero-crossing Doppler shifts, the asymmetry of observed V profiles is 
used for the diagnostics of motions in fiux tubes. This asymmetry detected for the first time 
in [37] from FTS observations near the disk center was "blue", i.e., the amplitude and area of 
the shortwave peak were greater, on the average, than in the longwave peak, the amplitude 
asymmetry being several fold greater than the area asymmetry. The general trends found in 
[57] have been confirmed by observations made with various resolutions, but the asymmetry 
magnitude slightly decreased with growing resolution. Numerous investigations demonstrated 
that the V profiles owe their asymmetry to the combined effect of the magnetic field and 
velocity field gradients inside fiux tubes and around them p9]. A large number of profiles 
with extremely large asymmetries were recently found in observations with resolutions below 
700 km. For example, three percent of all Fe I A 630.25 nm line profiles observed in the network 
regions and in faculae have such asymmetries [36j. Atypical V profiles, especially those which 
have one wing only, are of special interest. They were shown in MHD simulations to appear 
mainly at the periphery of magnetic elements [26]. Extremely large asymmetry can be caused 
by strong downfiows below the canopy near fiux tube boundaries. 

So, polarimetric observations with high spatial resolutions demonstrated that the dynamical 
characteristics of magnetic elements have a much wider range than that assumed before. What 
can we expect from observations in which the structure of magnetic elements will be resolved? 
We tackled this question with the use of the results of the 2D MHD simulation of magnetic 
elements and the synthesis of the Stokes profiles of the Fe I A 630.25 nm line. 

3 Two-dimensional MHD models 

The sequence of MHD models we use here was described in detail in [6], [71 [9]. The simulated 
region 3920x1820 km in size had 112 vertical columns (rays) with a 35-km spatial step. The 
starting magnetic field had bipolar decreasing with growing height, the mean unsigned field 
strength was 5.4 mT. The magnetic field was evolving in the course of two-hour simulation, 
and its unsigned strength grew to 50 mT on the average. Strong fiux tubes are formed after a 
lapse of 50 min of simulation. They disintegrate and form again in the simulation region. In 
one of the snapshots displayed in Fig. 1 two kilogauss fiux tubes of different polarities can be 
seen to approach each other and to destroy one the other a short time later. The horizontal 
size of flux tubes at the level of visible surface (logr = 0) varies from 35 km to 350 km in 
the course of evolution, and the maximum strength of their flelds varies correspondingly from 
40 mT to 250 mT. The largest area (about 50 percent) is occupied throughout the simulation 
time by flux tubes with diameters of 80-180 km and fleld strengths of 100-200 mT [6j. 

Our study is based on a sequence of the snapshots of half an hour's duration after 94 
minutes of simulation with a 30-s step. This is a period when strong flux tubes are intensely 
evolving. Two flux tubes of opposite polarities are often formed in the simulation region. 
As a rule, flux tubes are accompanied by nearby weak flelds of opposite polarity. In upper 
lateral layers of flux tubes hot regions appear as a result of deceleration of strong granulation 
flows near the tube walls. These hot regions may be analogs of observed bright points [1]. 




Figure 1: A snapshot from 2D MHD simulation jGj at moment 95.5 min: a) velocity field (arrow 
length is proportional to velocity magnitude) and isotherms (from top to botton) corresponding 
to 4000 K (thin line), 5000 K (thick line), 6500 K (dotted line), and 10000 K (thin line); b) field 
lines and magnetic field strength (shading density is proportional to field strengths of 1, 20, 
50, 100, 150 mT); c) field fines and field polarity (grey area shows positive polarity and white 
area negative polarity). Vertical axis) geometric height if, horizontal axis) column number N 
in the MHD model; distance between the columns is 35 km. 
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Figure 2: Ratio of V" profiles of Fe I lines 
AA 525.0 and 524.7 nm as a function of 
the distance from V^-profile center. Solid 
curves were obtained from MHD mod- 
els: 1) for region with one subkilogauss 
flux tube, 2) for region with two kilogauss 
flux tubes, 3) for axis of kilogauss flux 
tube; plusses) observations ^28j outside ac- 
tive regions; curve 4) observations [IH] for 
strong faculae; curve 5) observations [l6] 
for weak faculae. 



Predominantly horizontal magnetic flelds are observed in the central parts of granular cells. All 
these peculiarities of the simulated magnetogranulation can be seen in the vertical section of 
the simulation region (Fig. 1). We assume that the selected snapshots sequence represents the 
quiet Sun regions with some network elements because the unsigned fleld strength distribution 
mode is 35 mT in the simulation region and the flux density is 0.2 mT 



4 Testing the MHD models 

The MHD models we use were tested in ^32j. Here we continue to test them by the method 
of fleld strength ratio measured in two spectral lines |44j. This well-known method was widely 
used for the diagnostics of strong magnetic flelds outside active regions on the Sun in low spatial 
resolution spectra (1-4") [IHl [161 [T71 HO] , ^-iid now it is applied to high-resolution spectra (< 1") 
[T7t [52] . The effectivity of the line-ratio method depends on the model atmosphere chosen for 
the fleld strength calibration, so that the method also needs testing with models which are 
more realistic than the two-component models commonly used for calibration pJTj. Applying 
this method to MHD models, we can test the models and the method itself. 

As in observations, we chose two Fe I lines - AA 524.71 and 525.02 nm. Their atomic 
parameters are nearly the same, and only the Lande factors differ (their ratio is 1.5). Therefore, 
the ratios of calibrated magnetograph signals in two lines - -B525.0/-B524.7 or ^525.0/(1-5^524.7) 
- should be insensitive to all atmospheric parameters except fleld strength, magnetic fleld 
inclination, and velocity fleld. As shown in |2Sl 113 SS], the "magnetic ratio" also depends 
on the position of the selected V-proflle section with respect to the proflle center - as the 
magnetograph slit approaches the line center, the ratio decreases. The line ratio also depends 
on the flUing factor in the resolved area and on lateral magnetic fleld proflle. Thus, the magnetic 
fleld strength found by the line-ratio method depends on the spatial and spectral resolutions 
of observed spectra. 

We determined the ratio of V signals as a function of distance to line center MLR{AX) = 
V525.o(^A)/[1.5V524.7(AA)] for the proflles in different simulation regions. The Doppler shifts 
of V proflles with respect to the central wavelength of / proflles were compensated. Figure 2 
illustrates three MLR{AX) relations. One relation was derived from proflles averaged over a 3" 
region with a relatively weak subkilogauss flux tube {Bq < 90 mT), the second relation refers to 



an area of the same size with two strong flux tubes with Bq < 190 niT, and the third relation 
is based on the profiles corresponding to the axis of a strong flux tube with Bq > 180 mT 
at the level logrs = (without spatial averaging). We also plotted the relations for the 
ratios -B525.0/-B524.7 (observations from [28]) and V525.o(AA)/[1.5V524.7(AA)] (FTS observations 
of regions with strong and weak faculae |1Z]). The magnetic ratio is greater at greater distances 
from line center and is greater than unity in distant line wings. At lower resolutions, the relation 
is less steep and the magnetic ratio is smaller (cf curves 2 and 3). The run of relation 2 (the 
region with two strong flux tubes and with low resolution) is in good agreement with observation 
data, and this supports our earlier inference [32j that the theoretical MHD models ||6j represent 
quite well the actual magnetic fields on the Sun. 

The run of the MLR{AX) relations derived for different magnetic fluxes from MHD models 
is the same as the run derived earlier with the use of two-component models [281 SS] • These 
relations roughly represent various cases of magnetic saturation and various long-wave shifts 
of V signals in different sections of the profiles of two lines. Kilogauss fields not only diminish 
the ^-profile amplitude but substantially change the profile shape as well. In strong fields the 
V525.0 profile is lower and wider than the V524.7 profile, and the ratio of magnetic saturations in 
two lines becomes reverse at a distance of 5-7 pm from the line centers: MLR becomes greater 
than unity. So, the MLR{AX) relations we derived clearly demonstrate that the magnetic 
fields measured by the line-ratio method heavily depend on the line section (AA) chosen for 
the measurements. 

The ratio of V^-profile amplitudes MLR = ay,525.o/(l-5av',524.7), which is independent of 
AA, is sometimes used in polarimetric observations [16]. We determined the amplitude ratio 
for 5824 V profiles calculated in every MHD model column (Fig. 3a). There is a correlation 
between the "magnetic ratio" and the amplitude av,525.o'- the ratio decreases, on the average, 
with increasing l^-profile amplitude, and this suggest that the magnetic field increases. We also 
determined MLR for profiles averaged with a 1" resolution. This spatial averaging markedly 
affected the l^-profile amplitudes and consequently the magnetic ratio (Fig. 3b ); the MLR 
scatter increased in every interval and the correlation between MLR and profile amplitude 
deteriorated. The number of atypically shaped profiles (with one wing or with several wings) 
was found to grow among weaker V profiles, and this can affect the accuracy with which mean 
V^-profile amplitudes are calculated. It should be noted that the Fe I lines A A 524.71 and 
525.02 nm have very narrow profiles which are highly sensitive to temperature, so that the 
effect of extreme asymmetry of weak V profiles is stronger for these lines, and this accounts for 
the greater scatter in MLR for small ^-profile amplitudes. 

We calculated the MLR values in order to transform the magnetic ratios into field strengths. 
This was done with the data for the MHD model shown in Fig. 1. We calculated MLR and 
effective heights of l^-profile peak formation for each model column. Having found the level 
at which the method locates the magnetic field, we determined the field strength for this level 
directly from the model and plotted the MLR-B relation (Fig. 3c). We found from this relation 
that a decrease of MLR from 0.94 to 0.82 (Fig. 3a) corresponds to a growth of B from 30 to 
62 mT (Fig. 3c). This means that, according to the line-ratio method, the field strength in the 
simulation region can be 60 ± 30 mT at the most, i.e., the method did not fix any strong fields 
about 150 niT appropriate for the MHD models we use. The averaged magnetic line ratio is 
0.88 ± 0.10, and the corresponding mean field strength in the simulation region is 43 mT, in 
good agreement with the simulation results for our MHD model sequence (40-50 mT) and with 
the results obtained in [52] by the line-ratio method applied to the V profiles of the same two 
lines observed in quiet regions with high spatial resolution. No significant deviations of MLR 
from unity were found in f52] , and the intra- network fields were concluded to have no kilogauss 
structures; the typical field strength was estimated at 20-50 mT. The coincidence of the data 
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Figure 3: Magnetic ratio MLR = 
ai^,525.o/(l-5ay,524.7): a) without averaging, 

b) with spatial averaging over 700 km. 
Standard deviations are shown for some 
intervals with the same number of points; 

c) dependence of MLR on magnetic field 
strength B derived from the data of the 
snapshot shown in Fig. 1. 



from ^52j with the mean field strength in the MHD models we used allows us to suggest that 
our simulation region is similar to quiet solar regions. 

Thus, the testing of MHD models by the line-ratio method with the use of the Fe I lines 
AA 524.71 and 525.02 nm showed a satisfactory agreement between the calculated and observed 
magnetic ratios as functions of the distance from the V^-profile center. This suggests that the 
MHD models ^ describe adequately the fine structure of actual magnetic fields of the quiet 
Sun. The testing also demonstrated that fields stronger than 100 mT with filling factors of 
1-5 percent cannot be detected by the line-ratio method in quiet regions on the Sun, since this 
method depends on spatial resolution. Nevertheless, the method can give a reliable estimate 
for the mean field strength in quiet regions. 



5 Results of analysis of synthesized V profiles 

The Stokes profiles of the Fe I line A 630.25 nm were calculated for every column in the MHD 
models by integrating the Unno-Rachkovskii equations for the polarized radiation transfer in the 
LTE approximation [21]. We chose this moderately strong photospheric line for the ^-profile 
analysis because it is most often used, together with the Fe I line A 630.15 nm, in polarimetric 
measurements [H |20l [221 |23l [36] . A telluric O2 line which is observed in the /-profile wing does 
not affect the V profile. The line A 630.25 nm is preferred to the magnetic Fe I lines AA 524.7 
and 520.5 nm not only because of its unblended V profile but also due to its high sensitivity 
to magnetic field and much smaller sensitivity to temperature irregularities. The V profile of 
Fe I A 630.25 nm is formed rather deep in the photosphere, at the level logrs = —1, on the 
average, and thus it can be calculated in the LTE approximation. It should be stressed that 
the dynamical characteristics of magnetic elements studied with the use of this line refer to the 
same photospheric level. This is particularly true for the velocity and magnetic field gradients 
which can suddenly change in highly inhomogeneous media, for example, at the periphery of 
compact magnetic features. That is why the asymmetry magnitude and sign, which attest 
to the existence of gradients, as well as the relationship between the asymmetry and other 
parameters found in our study can differ from the results obtained with the use of other lines. 
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We compared our calculations to the FTS observation data kindly made available by J. Sten- 
flo and S. Solanki. These observations were made in quiet network elements and in active facula 
regions at the McMath telescope in 1979; the Fourier spectrograph had a high spectral reso- 
lution (420 000) and low spatial and time resolutions (10" and 35-52 min, respectively) ^7\. 
We used only the V profiles of Fe I and Fe II lines in the wavelength range from 445.0 to 
557.0 nm. The observed V profile parameters were calculated by Solanki's codes, and the 
absolute zero-crossing shifts were determined by the method proposed in [5]. 

We analyzed the following \^-profile parameters. 1. Mean amplitude of blue (b) and red (r) 
wings ay = (|ab| + \ar\)/2. 2. Doppler zero-crossing shift with respect to the absolute wave- 
length (transformed into the line-of-sight velocity Vz by the standard formula). 3. Amplitude 
asymmetry 6a = {\ab\ — |ar|)/(|ob| + |flr|) and area asymmetry 6A = \Ab\ — \Ar\)/{\Ab\ + \Ar\). 
The blue and red wing amplitudes (a;, and a^) and the zero-crossing position were calculated by 
fitting the corresponding profile sections to polynomials. The areas Af, and A,, were obtained by 
integrating the V/Ic profile from zero-crossing to a 0.5-percent level in the red and blue wings. 

By analogy with the observation data from [36] , we used only the profiles with amplitudes 
greater than 0.15 percent and regularly shaped profiles, i.e., only those with two wings of 
different signs and with one zero-crossing. The number of such profiles was 3755; the rest 
were either very weak or abnormal in shape (symmetric, with one wing only, or with several 
zero-crossings), and they were excluded from the analysis. 

Mean V-profile amplitudes. Figure 4 (upper part) demonstrates the ay distribution 
obtained from V^-profiles calculated without spatial averaging and with a 1" averaging. When 
the histogram for all profiles is compared to the histogram for regularly shaped profiles (without 
strong anomalies), one can see that weak profiles are less numerous in the second case due 
to highly asymmetry of weak profiles. Spatial averaging of profiles essentially changes their 
distribution. Recall that the maximum diameter of simulated flux tubes is about 350 km, so 
that the profiles averaged over a 700-km area do not resolve the tube structure, and profiles 
with amplitudes above 10 percent are absent in the histogram. 

In the analyses of observed Stokes profiles, ay is most often used as a substitute for filling 
factor or as an indicator of the longitudinal magnetic field. The statistical relations for ay 
shown in Fig. 4 (lower part) were averaged over the intervals of the field strength B and 
field inclination 7, they were obtained directly from the MHD models for the optical depth 
logTs = — 1 and from the continuum intensity contrast Ic/ < Ic >■ Here Ic is the continuum 
intensity calculated for a specific model column and < Jc > is the intensity averaged over 
the entire simulation region in the model. There is a close correlation between B and ay as 
well as between 7 and ay. The relationship between Ic/ < Ic > and ay is more complicated 
— the weakest \^-profiles correspond to high-contrast areas, the strongest profiles correspond 
to low-contrast areas, while moderately strong profiles represent areas with slightly changing 
contrast. 

So, the nearly linear ay-B relation found suggests that ay can be used as a magnetic field 
indicator in the analyses of Stokes profiles. 

Zero-crossing shifts. The line-of-sight velocities (Vz) of the magnetized plasmas fiows only 
are derived from the Doppler shifts of V profiles or the so-called zero-crossing shifts. Positive 
(red) shifts point to downfiows and negative (blue) shifts to upfiows. Figure 5 (upper part) 
displays the Doppler shifts converted into V^ as functions of three principal parameters derived 
usually from the observed ^-profiles: ay (indicator of magnetic field strength), magnetic vector 
inclination 7, and intensity contrast Ic/ < Ic >■ The inclination angle was derived from the 
amplitude ratio for Stokes profiles: tan^7 ^ Q"^ + U'^y^'^/V'^. 

The distribution of V^s (Fig. 5, left) points to a well-marked asymmetry between upfiows and 
downfiows - a fundamental property of magnetoconvection in a matter with frozen magnetic 
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Figure 4: Upper plots) distributions of V profile amplitudes ay from 52 MHD snapshots: 
at the left) for all 5824 profiles {oy = 7%); at the center) with abnormal profiles excluded 
(3755 profiles, ay = 9%); at the right) for profiles spatially averaged with a 700- km step (243 
profiles, ov = 4%). Lower plots) unsigned field strength, field inclination angle, and contrast 
of continuum intensity as functions of amplitude ay- 
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Figure 5: Vz histograms derived from the zero-crossing shifts of V profiles, and the V^ depen- 
dence on ay, 7, and Id < Ic >■ upper plots) all profiles, lower plots) profiles averaged with a 
700-km step. Dotted line and squares) FTS observation data. 
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Figure 6: Oscillations derived 
from 2D MHD model of solar 
granulation, a) Vertical velocity 
Vz (curve 1), gas pressure fluc- 
tuations 5P (curve 2). Dashed 
straight line) mean velocity of 
0.6 ± 1.1 km s-^ b) V, (curve 1), 
unsigned magnetic field strength 
5B at the level logrs = (curve 
3). Vz, SP, and 6B are averaged 
over whole simulation region. 



field. The velocity distribution maximum lies at —1 km s~^ approximately. A mean velocity 
of 0.53 km s~^ and a wide velocity range from —3 km s~^ to 9 km s~^ suggest that downfiows 
dominate in the magnetized plasma. There is a correlation between velocity and principal 
parameters of V profiles . The downfiow (positive) velocity increases to 3 km s^^, on the 
average, with increasing ay (-B), with decreasing 7 and Ic/ < Ic >, and this suggests that the 
most rapidly moving plasma resides in strong (120 mT) nearly vertical (10°) flux tubes with 
low contrast (0.7), which corresponds to darker intergranular lanes. The above correlations are 
in accord with the physical properties of convective motions which affects the distribution of 
LOS velocities and magnetic fields in photosphere (see velocity field of the convective motions 
in the snapshot. Fig. 1). 

The spatial averaging of profiles decreases the fraction of very weak profiles with negative 
velocities and very strong profiles with high positive velocities (the lower part of Fig. 5). The 
spatial averaging increased the mean value of Vz to 0.72 km s~^. All relations in the lower part 
of Fig. 5 preserved the trends pronounced in the upper plots. 

Our results agree satisfactorily with the observations of magnetic regions outside activity 
centers analyzed in [36], where the mean velocity was estimated at 0.73 km s~^ and the velocities 
ranged from —6 to 6 km s~^. In [19] Vz was 0.25 km s~^ and the range was from —3 to 5 km s~^. 
We also plotted the FTS observation data from |17] [Vz = —0.04 km s~^) in Fig. 5. Our results 
are in good agreement with observed velocities in the regions of small ay- 

Oscillations. Although wave processes in spatially unresolved flux tubes are difficult to 
detect directly, oscillatory motions with a 5-minute period along flux tubes have been reported 
[38]. There is also some evidence for oscillations with a shorter period [51]. Observations with 
high spatial resolutions [1] revealed velocity oscillations with an amplitude of about 1 km s~^. 
Theory predicts a wide variety of wave modes in solar flux tubes and their importance in the 
heating of the chromosphere and corona. According to the 2D MHD simulation [B] velocities 
of vertical flows in magnetized plasmas can vary in a wide range (see Fig. 5). It points to an 
oscillatory instability of magnetogranulation regions. Here we examined the temporal changes 
of velocity Vz derived from shifts of V profiles averaged over in the whole simulation region 
for the complete two-hour sequence of MHD models with a 1-minute interval. Fig. 6a shows 
obtained vertical velocities as functions of time. We also calculated the power spectrum of 
velocity oscillations (Fig. 7a), where two power peaks can be clearly seen in the bands at 5 min 
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(2 < z/ < 4.5 mHz) and 3 min (4.5 < u < 7 niHz). The power maximum lies at i/ ~ 2.8 mHz 
(5.9 min), and a smaller peak lies at z/ ~ 4.7 mHz (3.5 min). A 99-percent confidence level 
shown by a horizontal line in the figure was calculated in accordance with [H] from the mean 
power in the band 5.5-8.3 mHz. The Nyquist frequency is equal to 8.3 mHz. 

Velocity oscillations with frequencies of 2.57, 3.88, and 5.58 mHz were also obtained in [12] 
from a 3-D simulation of the solar surface convection in the absence of magnetic field. The 
p-mode oscillations were assumed to be excited by random nonadiabatic pressure fluctuations 
near the Sun's surface. Such fluctuations are produced by the radiative cooling which can 
locally deviate for a short time from equilibrium with the heating produced by convective 
motions. Figure 6a shows the gas pressure fluctuations 6Pg calculated for our models. Here 
Pg is the pressure averaged over one model at a specific point in time and \Pgi is the pressure 
averaged over all models in a time interval of 120 min. Comparison of the velocity and pressure 
oscillations suggests that there should be a relationship between them. In the power spectrum 
calculated for pressure oscillations (Fig. 7b) there are power peaks at the same frequencies as 
in the spectrum of velocity oscillations in the 5-min and 3-min bands. 

To examine the influence of magnetic field on the oscillations of vertical velocities, we 
obtained the field strength fluctuations 6B (Fig. 6b) from the simulation data and calculated 
their power spectrum (Fig. 7c). First of all we point out a constant increase of the magnetic 
field and its well-marked oscillations associated with the evolution and disruption of strong 
flux tubes in the simulation region. A strong peak in the field strength oscillation power was 
found at 1/ = 0.83 mHz (20 min) and a smaller peak was found at 1.25 mHz (13 min). When 
comparing the magnetic field oscillations with vertical velocity oscillations, we found that the 
negative velocity component almost disappeared and the positive component increased when 
the magnetic field was stronger, that is, the upflows were slower and downflows faster. This 
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Figure 8: Distributions of amplitude asymmetry 5a and area asymmetry SA of V profiles, 
scatter of asymmetry parameters and their correlations: upper plots) all profiles, lower plots) 
profiles averaged with a 700-km step. Dotted line and square) FTS observation data. 



means that the velocity oscillations are nonlinear. Such quasi-oscillatory motions in flux tubes 
were predicted in [12]. Nonlinear oscillations were considered there to be a superposition of 
upfiows and downfiows with different amplitudes and durations. The upflow velocities were 
estimated there at —0.5 km s~^ and downfiow velocities at about 2 km s~^. Our results for 
the moments of magnetic field strengthening are very close to these data. As seen in Fig. 6a, 
the magnetic field modulates the velocity oscillations. Based on the results of MHD simulation 
[6], we attribute this effect to convective instability in intergranular lanes. About 10-12 min 
after the formation of a strong flux tube begins, conditions favorable for convective collapse set 
up, and the collapse produces a kilogauss field. Then the flux tube is in a quasi-stable state 
for some time (< 10 min), it is narrower and closer to the vertical, its temperature grows, the 
field is stronger, and the Wilson depression correlates with the vertical velocity. Oscillatory 
downward motions become more intense in the flux tube region. After each abrupt downward 
shift of magnetic configuration, there comes a period of a slacker convective collapse. These 
oscillations persist until the strong evacuation in the upper part of the flux tube begins to 
destroy the tube (i.e., until the beginning of reverse convective collapse). As this takes place, 
downfiows are replaced by upfiows with supersonic velocities. The flux tube rapidly dissipates 
in the course of several minutes. The size and the inclination of the flux tube increase, the 
magnetic field and the temperature diminish, the Wilson depression also diminishes. As the 
dynamical condition changes, the inclination of flux tubes also changes with time — it oscillates. 
This can give rise to waves propagating along the flux tubes |43j . 

So, the convective and superconvective instability processes going in strong flux tubes bring 
about a nonlinearity in the oscillatory motions of magnetized plasmas in the simulation region. 

Asymmetry. Figure 8 demonstrates the distribution of the amplitude and area asymme- 
tries, 6a and 6A, and their correlations. For unaveraged profiles, the amplitude asymmetry 
varies in a wide range (±70 percent) at small amplitudes. The area asymmetry has a similar 
distribution, but the scatter is smaller. The amplitude asymmetry correlates with the area 
asymmetry, at positive values they vary in nearly direct proportion. This means that both 
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Figure 9: Amplitude asymmetry as a function of mean amplitude of V profile, vertical velocity, 
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averaged with a 700-km step. 



asymmetries are caused by the same factors - the velocity field and magnetic field gradients. 
Mean asymmetries are close to zero {5a = —1, 6 A = 1 percent for unaveraged profiles and 
6a = 3, 6 A = —2 percent for averaged ones). Spatial averaging drastically weakens the correla- 
tion between 6a and 6 A and changes the inclination of the relation — the amplitude asymmetry 
becomes more positive. This suggests that the spatial averaging is also a cause of observed pro- 
file asymmetry, the amplitude asymmetry being more sensitive to the averaging. This may be 
the reason why the amplitude asymmetry measured from observations is greater, as a rule, than 
the area asymmetry: 6a = 15 percent and 6 A = 6 and 4 percent for the Fe I lines AA 630.15 and 
630.25 nm [361 [TT]. Figure 8 also shows the asymmetries found in FTS observations |17], they 
are much greater, on the average, than those obtained in our calculations. This discrepancy 
seems to be caused by a substantial difference in the spatial and temporal averaging and by 
different activity levels in the observed and simulated regions. 

The observed area asymmetry is a factor of 2-3 greater than the amplitude asymmetry. 
This difference is yet to be explained. Our results suggest that it is caused by insufficient 
spatial resolution of observations and by atmospheric effects. When the quality of observations 
and the spatial resolution are upgraded, the relative number of strong V profiles observed in 
the kilogauss flux tube will increase. Our calculations demonstrate that the asymmetry in 
such profiles is much smaller because there are no sudden velocity and field strength gradients 
inside flux tubes as compared to their periphery or the regions of weak turbulence fields of 
mixed polarity. These properties of strong V profiles also account for a small scatter in their 
asymmetries, while the scatter for weak V profiles can be as large as 70 percent. 

We analyzed the relationship between the ^-profile asymmetry and the amplitude ay, ve- 
locity Vz, inclination angle 7, and contrast (Fig. 9). The correlation of 6a with these parameters 
is rather weak; the only feature worthy of notice is an appreciable growth of positive asymmetry 
with the velocity of upfiows while downflows are more often associated with negative amplitude 
asymmetries. 

Center-to-limb variation. Investigations of the center-to-limb variations in the V^-profile 
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shifts and asymmetries involve a huge amount of calculations, and so we selected ten snapshots 
only. Figure 10 shows the velocity Vz as functions of their position on the disk (/i = cos 6). Vz 
is positive (predominance of downflows) outside the disk center, it amounts up to 2 km s~^, on 
the average, at /x = 0.9 and changes almost not at all to the limb. 

Spatial averaging slightly affects the Vz-fi dependence. According to observations ^^, 
downflows dominate at the disk center, and the profile shifts are close to zero at the limb. Our 
results can differ from observations for various reasons. We suppose that the velocity increase 
obtained in our calculations is a result of the canopy effect in the upper parts of flux tubes, 
where the profiles calculated for cos 6^1 are formed. According to the 2D MHD simulations 
used in this paper, the strong flows of matter in the upper layers are predominantly inclined, 
they are directed slightly downwards, to the centers of intergranular lanes (see Fig. 1). Another 
reason for the discrepancy can be related to some peculiarities of the 2D MHD models - the 
upper boundary conditions imposed on the upper layers in the models do not correspond to 
actual conditions in the upper photosphere. This may affect the calculations of the line profiles 
far from the disk center, these profiles being formed higher in the photosphere than those which 
are located closer to the disk center. At the same time the accuracy of observations deteriorates 
closer to the limb. It should also be noted that there is observational evidence of high horizontal 
velocities (up to 2 km s~^) [51], and they are in better accord with our results. 

The center-to-limb variations of mean asymmetry 6a, particularly the asymmetry for aver- 
aged profiles, agree satisfactorily with observation data. In FTS observations fJS] the amplitude 
asymmetry diminished from 10-15 percent at the disk center to zero at the limb. In [TT] the 
mean asymmetry 6a varied from 15 percent (center) to —5 percent (limb), with the intersection 
point near cos^ = 0.7. Similar results were obtained in [23], but 6a and 6 A slightly increased 
at cos^ = 0.86 and then monotonically diminished to the limb. 

The area asymmetry in our calculations is small and predominantly negative (Fig. 10), while 
in the FTS observations in ^3| the asymmetry is positive at the center, diminishes to the limb 
and becomes negative. Negative area asymmetries were observed in [21] at the disk center in 
regions with small filling factor. 
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One can see in Fig. 10 that the spatial averaging of profiles affected the most the amplitude 
asymmetry, especially at the disk center. We wish to stress once more that the amplitude 
asymmetry is the parameter most sensitive to spatial averaging. 

6 Conclusion 

We used the results of two-dimensional MHD simulation of solar granulation to investigate the 
motions of matter in magnetic elements with high spatial resolution (35 km). The analysis of 
the V profiles of the Fe I A 630.25 nm line synthesized in snapshots of the MHD simulations 
allowed us to make some inferences as to the vertical velocity of magnetized plasma motions 
in photospheric regions with a mean magnetic flux density of 0.2 mT and mean unsigned field 
strength of 35 mT. 

1. The mean velocity is 0.5 ± 2 km s~^ with a scatter from —3 to 9 km s~^. Downfiows 
with a mean velocity of 3 ± 2 km s~^ dominate in intergranular regions inside strong vertical 
flux tubes, while upflows with a mean velocity of 0.5 ± 2 km s~^ are typical of granular regions 
outside flux tubes. 

2. There is a noticeable correlation between the velocity and the amplitude and area asym- 
metries of V profiles. The positive asymmetry increases, on the average, with upfiow velocity. 
The negative asymmetry occurs most often in the profiles formed at the sites with dominant 
downfiows. The mean asymmetry of amplitudes and areas is about 1 percent with a scatter of 
70 percent in weak V profiles and 10 percent in strong ones. 

3. The mean velocity fiuctuates nonlinearly in time. There are two peaks in the velocity 
power spectrum — a stronger peak in the 5-min band and a much weaker one in the 3-min 
band. Velocity fiuctuations correlate with gas pressure fiuctuations at the surface level, and they 
also depend on magnetic field fiuctuations. The power spectrum of field strength fiuctuations 
has two peaks at periods of 20 and 13 min. Periodic variations of magnetic field strength 
are related to intensification and dissipation of magnetic fields in strong fiux tubes. When 
the field strength increases, oscillatory downward motions are intensified, and this results in a 
nonlinearity of mean velocity oscillations. 

4. The mean velocity substantially changes when going from the solar disk center to the 
limb. These changes occur due to intense slopping motions in the upper layers of simulation 
region. The mean velocity increases to 2 km s~^ at a small distance from the disk center 
(cos^ = 0.9), and then it slightly changes to the limb. The amplitude asymmetry decreases 
and changes its sign from positive to negative when going from the center to the limb of the 
solar disk. The negative area asymmetry is greater at the limb. 

5. The mean velocity depends on spatial averaging of profiles. At lower resolutions the 
^-profile amplitudes and other parameters substantially change, and they are closer to the 
observed parameters. The asymmetry of V profiles is the parameter most sensitive to spatial 
averaging, the amplitude asymmetry growing much faster at lower resolutions than the area 
asymmetry. This may be a reason why observed amplitude asymmetries are greater than area 
asymmetries. 
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